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The  shunt  current  loss  is  one  of  main  factors  to  affect  the  performance  of  the  vanadium  redox  flow  battery, 
which  will  shorten  the  cycle  life  and  decrease  the  energy  transfer  efficiency.  In  this  paper,  a  stack-level 
model  based  on  the  circuit  analog  method  is  proposed  to  research  the  shunt  current  loss  of  the  vanadium 
redox  flow  battery,  in  which  the  SOC  (state  of  charge)  of  electrolyte  is  introduced.  The  distribution  of  shunt 
current  is  described  in  detail.  The  sensitive  analysis  of  shunt  current  is  reported.  The  shunt  current  loss 
in  charge/discharge  cycle  is  predicted  with  the  given  experimental  data.  The  effect  of  charge/discharge 
pattern  on  the  shunt  current  loss  is  studied.  The  result  shows  that  the  reduction  of  the  number  of  single 
cells  in  series,  the  decrease  of  the  resistances  of  manifold  and  channel  and  the  increase  of  the  power  of 
single  cell  will  be  the  further  development  for  the  VRFB  stack. 
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1.  Introduction 

The  vanadium  redox  flow  battery  (VRFB)  is  one  of  most  promis¬ 
ing  large-scale  storage  technologies  to  meet  the  requirement  of  grid 
load-smoothing  and  smooth  output  of  renewable  energy  sources, 
due  to  its  characters  like  high  performance,  long  cycle  life  and 
flexible  design  [1-3].  The  cell  stack  is  assembled  by  several  single 
cells,  which  are  connected  in  series  by  bipolar  plates  to  meet  the 
voltage  requirement  [4].  Each  single  cell  shares  the  anode/cathode 
electrolyte  with  the  anode/cathode  manifold.  The  single  cells  are 
connected  with  the  manifold  and  their  channels.  Because  of  the 
conductivity  of  electrolyte  and  non-zero  electrical  field  poten¬ 
tial  gradient,  the  shunt  current  will  be  exited.  Two  problems  are 
induced  by  the  presence  of  shunt  current,  one  is  the  corrosion  of 
the  key  materials  of  cell  stack,  which  will  shorten  the  cycle  life;  the 
other  is  the  capacity  loss,  which  will  decrease  the  energy  transfer 
efficiency  of  the  cell  stack  [4].  Thus,  more  efforts  should  be  done  to 
reduce  and  even  eliminate  the  shunt  current  loss. 

Up  to  now,  most  researches  about  VRFB  were  focused  on  the 
key  materials  like  electrode,  current  collector  and  ion  exchange 
membrane.  Very  few  reports  are  related  to  the  shunt  current  of  the 
vanadium  redox  flow  battery.  White  et  al.  proposed  a  method  based 
on  the  circuit  analog  model  for  predicting  shunt  current  in  stack  of 
divided  plate  cells  [5].  Schaeffer  et  al.  presented  a  stack-level  model 
based  on  the  electrical  circuit  to  calculate  the  shunt  current  of  fuel 
cells  [6].  Henquin  developed  a  simplified  mathematical  model  to 
calculate  the  current  distributions  in  bipolar  electrochemical  reac¬ 
tors  [7].  In  the  available  references,  the  sensitive  analysis  of  shunt 
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current  and  the  shunt  current  loss  are  not  considered.  Moreover, 
the  literatures  are  lack  of  the  research  on  the  shunt  current  of  VRFB. 

In  this  paper,  a  stack-level  model  based  on  the  circuit  analog 
method  is  proposed  to  investigate  the  shunt  current  of  VRFB,  in 
which  the  SOC  (state  of  charge)  of  electrolyte  is  introduced.  The 
model  is  validated  by  the  given  experimental  data.  The  distribu¬ 
tion  of  shunt  current  is  described  in  detail.  The  sensitive  parameters 
affecting  shunt  current  such  as  voltage  of  the  single  cell,  the  number 
of  single  cells  in  series  and  the  resistances  of  manifold  and  chan¬ 
nel  are  considered.  The  shunt  current  loss  in  the  charge/discharge 
cycle  of  cell  stack  is  predicted  with  the  given  experimental  data. 
The  effect  of  charge/discharge  pattern  on  the  shunt  current  loss  is 
studied.  The  further  development  for  the  VRFB  is  suggested. 


2.  Model  description 

Shunt  current  of  the  vanadium  redox  flow  battery  can  be  simpli¬ 
fied  as  shown  in  Fig.  1  based  on  the  circuit  analog  model,  since  the 
electrolyte  inlet  loop  and  outlet  loop  of  the  cell  stack  is  symmetrical. 

Where  Re  is  the  resistance  of  single  cell,  RA  and  Rc  is  the  chan¬ 
nel  resistance  for  anodic  and  cathodic  electrolyte,  respectively;  RMA 
and  RMc  is  the  manifold  resistance  for  anodic  and  cathodic  elec¬ 
trolyte,  respectively;  V0  is  the  open-circuit  voltage  for  each  cell; 
Vn  is  the  voltage  for  the  nth  single  cell;  k  and  I<  is  the  channel  and 
manifold  current  for  anodic  electrolyte,  respectively;  l  and  L  is  the 
channel  and  manifold  current  for  cathodic  electrolyte,  respectively; 
in  is  the  current  passed  through  the  nth  single  cell;  It  is  applied 
current. 
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Fig.  1.  The  circuit  analog  model  for  a  cell  stack  of  N  single  cells. 


Application  of  Kirchhoffs  law  to  this  symmetric  network  at  the 
nth  cell  analog  results  in  a  set  of  five  linear  constant  coefficient 


algebraic  equations. 

Kfi+l  —  Kn  —  l<n+ 1  =  0 

(i) 

Ln+ 1  ~  Ln  —  ln+ 1  =  0 

(2) 

*n+l  —  in  +  2/<n+i  +  2/n+i  =  0 

(3) 

t^MC^n  +  Rc(in  —  ln+ 1 )  = 

(4) 

RmaKu  +  RaU<u  —  kn+1  )  —  Vn 

(5) 

For  the  first  cell  and  the  Nth  cell  in  cell  stack,  there  are  special 
relationships  as  followed: 


/<!=/<!  (6) 

lN  =  -Ln  (7) 

h  =k-  2/<i  (8) 

Resistances  of  channel  and  manifold  for  anolyte  and  catholyte, 
respectively,  can  be  calculated  by  [9]: 

Ri  =  Plf.  O) 


V?e//,  denotes  the  standard  cell  potential  at  50%  SOC.  The  value  of  it  is 
equal  to  1 .4  V  (25  °C)  which  is  derived  from  Haddadi’s  experimental 
result  [9].  The  relationship  between  charge  (or  discharge)  time  t  and 
SOC  is: 

SOC  =  So  +  -t(Se-S0)  (12) 

to 

where  t0  is  the  total  charge  (or  discharge)  time.  S0  and  Se  repre¬ 
sent  the  values  of  SOC  at  the  beginning  and  the  end  of  charge  (or 
discharge).  The  over  potential  of  each  cell  depends  on  the  inner- 
resistance  Re,  and  current  passed  through  the  nth  cell  in: 

F]n  —  in  •  Re  (13) 

The  output  terminal  voltage  is  a  summation  of  voltage  of  each 
cell: 

N 

Vstack  =  (14) 

n=l 

The  shunt  current  loss  during  charge  and  discharge  process  can 
be  calculated  by: 

rk  N  rtd  N 

Qf=  ^(t)dt+/  (15) 

n= 1  n= 1 


p ,  denotes  electrolyte  resistivity.  represents  length  of  the  nth 
electrolyte  pipe,  A/  is  cross  sectional  area  of  the  nth  electrolyte  pipe. 
In  this  paper,  the  voltage  of  the  nth  cell  is  related  to  SOC  (charge 
of  states)  and  over-potential.  This  means  the  model  can  simulate 
the  variation  of  output  terminal  voltage  in  the  whole  charging  and 
discharging  process.  It  is  important  that  the  characteristic  of  shunt 
current  changing  with  time  can  be  well  captured.  The  voltage  of  the 
nth  cell  is  represented  as  [8]: 


Vn  =  V0  +  rjn 


(10) 


V0>  denotes  the  open-circuit  voltage,  which  is  a  function  of  SOC: 


Vo  =  V°„  +  §  in 


[SOC]2 
[1  -  SOC]2 


(11) 


where  Pn  (or  iff)  is  the  total  shunt  current  in  the  process  of  charge 
(or  discharge)  for  the  nth  cell,  which  is  a  summation  of  absolute 
values  of  shunt  currents  in  manifold  and  channel,  respectively.  tc 
(or  td )  is  the  charge  (or  discharge)  time. 

3.  Results  and  discussion 

3.1.  Experiment  and  results 

A  kilowatt-class  VRFB  stack  is  employed  to  investigate  the 
shunt  current.  The  stack  is  manufactured  by  10  VRFB  sin¬ 
gle  cells.  The  effective  area  of  each  electrode  is  476  cm2.  The 
charge/discharge  performance  of  the  VRFB  stack  is  conducted  by 
using  a  charge-discharge  controller  ( Arbin  Instruments  Corp.,  USA). 


Table  1 

The  experimental  results. 


Applied  current  density/mAcm-2 

Charge  capacity/Ah 

Discharge  capacity/Ah 

So/% 

Sel% 

Charge  time/h 

Discharge  time/h 

40 

766 

695 

0.9 

78.4 

4.04 

3.66 

60 

707 

659 

1.5 

71.8 

2.48 

2.31 

80 

607 

573 

2.9 

63.6 

1.6 

1.51 
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Fig.  2.  Comparison  of  experimental  and  simulated  charge/discharge  curves  for  dif¬ 
ferent  applied  current  densities. 


Fig.  3.  Distributions  of  shunt  current  in  manifold  and  channel  during  charge  and 
discharge  at  50%  SOC  for  applied  current  density  of  60  mAcm-2. 


Charge 


The  upper  limit  voltage  of  charge  is  15.5  V.  The  lower  limit  voltage 
of  discharge  is  1 0  V.  The  experimental  results  are  shown  in  Table  1 . 

3.2.  Validation  of  the  model 

The  resistances  of  manifold,  channel  and  single  cell  of  the 
kilowatt-class  VRFB  stack  are  1.23,  365.89  and  0.0036  £2,  respec¬ 
tively.  The  charge/discharge  curve  is  simulated,  which  is  compared 
with  the  experimental  curve.  As  shown  in  Fig.  2  the  two  curves  are 
well-fitted  with  each  other.  The  slight  discrepancies  between  them 
are  caused  by  inaccurate  evaluation  of  the  SOC  at  the  beginning  and 
the  end  of  charge  or  discharge. 

3.3.  Distribution  of  shunt  current 

In  Fig.  3,  the  maximum  value  of  shunt  current  in  manifold  occurs 
in  the  middle  of  stack,  and  decreases  towards  two  ends.  The  vari¬ 
ation  demonstrates  a  symmetrical  profile.  The  maximum  shunt 
current  in  channel  occurs  at  two  ends  of  cell  stack.  The  minimum 
occurs  in  the  middle,  which  is  nearly  equal  to  zero.  The  sum  of  shunt 
currents  in  channel  is  equal  to  zero. 

Fig.  4  shows  the  distributions  of  single  cell  voltage  and  internal 
current  during  charge  and  discharge  at  50%  SOC  for  applied  current 
density  of  60  mA  cm-2.  The  distribution  of  single  cell  voltage  shows 
the  “U”  profile  during  the  charge  and  discharge  process.  The  inter¬ 
nal  current  passed  through  single  cell  shows  the  “U”  profile  during 
the  charge  process,  however,  shows  the  opposite  profile  during  the 
discharge  process. 

3.4.  Sensitive  analysis  of  shunt  current 

The  shunt  current  is  dependent  on  the  voltage  of  the  single  cell, 
the  number  of  single  cells  in  series  and  the  resistances  of  mani¬ 
fold  and  channel.  The  shunt  current  during  charge/discharge  cycle 
for  applied  current  density  of  80 mAcm-2  is  shown  in  Fig.  5.  The 
variation  is  directly  proportional  to  the  variation  of  the  single  cell 
voltage. 

With  the  increase  of  the  number  of  single  cells  in  series,  the  cir¬ 
cuit  loops  generating  shunt  current  increase.  In  order  to  analyze 
the  influence  of  the  number  of  single  cells  in  series,  the  shunt  cur¬ 
rent  of  the  stack  assembled  by  20  single  cells  in  series  is  calculated. 
Compared  with  the  internal  current  of  the  stack  assembled  by  10 
single  cells  in  series  as  shown  in  Fig.  6,  the  internal  current  becomes 
more  than  thrice  at  the  same  SOC.  So  the  reduction  of  the  number 
of  single  cells  in  series  is  an  effective  method  to  reduce  the  shunt 
current. 


Discharge 


Fig.  4.  Distributions  of  single  cell  voltage  and  internal  current  during  charge  and  discharge  at  50%  SOC  for  applied  current  density  of  60  mAcm-2. 
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Fig.  5.  Shunt  current  during  charge/discharge  cycle  for  applied  current  density  of 
80  mA  cm-2. 


The  shunt  current  and  the  resistances  of  manifold  and  channel 
are  in  inverse  variation.  The  effect  of  resistance  on  the  shunt  current 
is  investigated  via  increasing  the  resistances  of  channel  and  man¬ 
ifold  by  50%,  respectively.  As  shown  in  Fig.  7  the  former  and  the 
later  shunt  currents  are  decreased  by  1.64%  and  32.57%,  respec¬ 
tively.  It  can  be  concluded  that  under  the  condition  of  increasing 
equal  percentage  of  resistance,  increasing  resistance  of  channel  is 
more  effective  to  decrease  the  shunt  current. 

3.5.  Shunt  current  loss 

The  shunt  current  loss  for  different  applied  current  densi¬ 
ties  in  the  given  experiment  is  calculated  as  shown  in  Fig.  8. 
With  the  increase  of  applied  current  density,  the  shunt  cur¬ 
rent  loss  and  the  percentage  of  the  shunt  current  loss  in 
the  total  coulomb  loss  both  decrease.  The  former  is  mainly 
because  of  the  reduction  of  the  charge/discharge  time,  and 
the  later  is  due  to  the  increase  of  the  power  of  the  stack. 
In  the  given  experiment  the  percentage  of  the  shunt  current 
loss  in  the  total  coulomb  loss  is  less  than  1 7%.  In  case  that  the  num¬ 
ber  of  single  cells  and  the  area  of  the  electrode  are  kept  constant, 
increasing  the  applied  current  density  is  an  effective  method  to 
decrease  the  shunt  current  loss. 


Fig.  6.  Comparison  of  internal  current  for  different  number  of  cells  in  series  at  50% 
SOC  for  applied  current  density  of  80  mAcrrr2. 
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Fig.  7.  Comparison  of  shunt  current  for  different  resistances  at  50%  SOC  for  applied 
current  density  of  80  mAcm-2.  For  tick  of  x-axis:  (a)  original;  (b)  resistance  of  man¬ 
ifold  increased  by  50%;  (c)  resistance  of  channel  increased  by  50%. 


Fig.  8.  Comparison  of  shunt  current  loss  for  different  applied  current  densities. 


As  described  above,  the  shunt  current  is  constant  at  the  given 
SOC  if  the  number  of  single  cells  in  series  and  the  resistances  of 
manifold  and  channel  are  specified.  The  higher  the  power  of  the 
single  cell,  the  lower  the  influence  of  the  shunt  current  loss  is.  So 
the  further  efforts  on  VRFB  stack  should  be  focused  on  reducing 
the  number  of  single  cells  in  series,  decreasing  the  resistances  of 
manifold  and  channel  and  increasing  the  power  of  single  cell. 


4.  Conclusion 

A  stack-level  model  based  on  the  circuit  analog  method  is  pro¬ 
posed  to  investigate  the  shunt  current  loss  of  the  vanadium  redox 
flow  battery,  in  which  the  SOC  (state  of  charge)  of  electrolyte  is 
introduced.  The  model  is  validated  against  the  stack  voltage  of  the 
given  experiment,  which  is  well-fitted  with  the  simulated  curve. 
The  sensitive  analysis  of  the  voltage  of  the  single  cell,  the  number 
of  single  cells  in  series  and  the  resistances  of  manifold  and  channel 
to  the  shunt  current  is  reported.  The  shunt  current  loss  of  cell  stack 
in  the  charge/discharge  cycle  is  calculated  with  the  given  experi¬ 
mental  data.  The  percentage  of  the  shunt  current  loss  in  the  total 
coulomb  loss  is  less  than  17%.  The  effect  of  charge/discharge  pat¬ 
tern  on  the  shunt  current  loss  shows  that  the  reduction  of  number 
of  single  cells  in  series,  the  decrease  of  the  resistances  of  manifold 
and  channel  and  the  increase  of  the  power  of  single  cell  are  the 
further  development  for  the  VRFB  stack. 
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The  VRFB  system  is  constructed  by  several  stacks  in  parallel,  in 
which  the  shunt  current  exits  at  the  same  time.  The  proposed  model 
can  be  extended  to  calculate  the  shunt  current  of  VRFB  system. 
The  results  and  validity  of  the  model  for  the  VRFB  system  will  be 
discussed  in  a  forthcoming  paper. 
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